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Abstract The random local mixture of a family of primitive substitution rules with 
noble mean inflation multiplier is investigated. This extends the random Fibonacci 
example that was introduced by Godreche and Luck in 1989. We discuss the struc- 
ture of the corresponding dynamical systems, and determine their entropy, an er- 
godic invariant measure and diffraction spectrum. 



1 Introduction 



Despite many open problems (including the famous Pisot substitution conjecture), 
the structure of systems with pure point diffraction is rather well understood ll5l [T3l . 
Due to recent progress HI |2l, also the situation for various systems with diffrac- 
tion spectra of mixed type has improved, in particular from a computational point 
of view. In particular, one can explicitly calculate the diffraction measure in closed 
form for certain classes of examples. Still, the understanding of spectra in the pres- 
ence of entropy is only at its beginning; compare |[Tl[4l and references therein. The 
purpose of this contribution is a further step into 'disordered territory', here via the 
analysis of mixed substitutions that are randomly applied at a local level. This is 
in contrast to global mixtures (which leads to 5-adic systems), which have no en- 
tropy. Local mixtures were investigated in [Sj, where the essential properties of the 
Fibonacci case were derived, along with first results on planar systems based on 
triangle inflation rules. Here, we extend the random Fibonacci system to the noble 
means family, and present the results from the point of view of dynamical systems. 
The entire family is still relatively simple because each individual member of a fixed 
noble mean family defines the same (deterministic) hull. Various generalisations are 
possible, but not discussed here. 
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2 Construction 

Let = {a,^?} be our two letter alphabet. For any fixed integer m > 1, we define a 
family of substitution rules by 



771. Z ■ 



a 



with < ; < m, and refer to each , as a noble means substitution (NMS). 
Each member of is a primitive substitution with Pisot inflation multiplier 
km = \ {m + Vm^ +4-) and algebraic conjugate = j{m— V^i^ +4). Each sub- 
stitution possesses a reflection symmetric and aperiodic two-sided discrete (or sym- 
bolic) hull X„j where the hull, as usual, is defined as the orbit closure of a fixed 
point in the local topology. Moreover, all elements of Jif^n are pairwise conjugate to 
each other which implies that, for fixed m, the hulls X,„ , are equal for < / < m. 

We now fix a probability vector {pQ,...,p^), that is p^ > and L'jLoP; = 1- We 
define the random substitution rule 

{C,„.o(fl) ^ ba"\ with probability 
: : where M ( 

Cm.mi^) = with probability /7„j, \^ ^. 

bi-^ a 

is its substitution matrix. We refer to i^„, as the random noble means substitution 
(RNMS). The application of occurs locally, which means that we decide sepa- 
rately on each letter a which of the m + 1 possible realisations we choose. In partic- 
ular, for each A: e N, Cml'^) is ^ random variable. As there is no direct analogue of 
a fixed point in the stochastic situation, we have to slightly adjust the notion of the 
two-sided discrete hull in this context. Note that aa is a legal word (see below for 
more) for all m, and consider 

•= {we si'^ I wis an accumulation point of (5^"^^,(fl|fl))^jgpj^}, 

where S denotes the shift. The two-sided discrete hull X,„ is defined as the smallest 
closed and shift-invariant subset of si'^ with X„, C X„,. It is immediate that X,„ is 
a superset of X^ Note that typical elements of X„, contain the sub word bb, which 
is absent in X^ The hull X„, is characterised by the property that it contains all 
elements of {a,^}^ that contain i^,„-legal subwords only (see below for more). 



3 Topological entropy 



In this section, we assume that all probabilities are strictly positive. We call a 
finite word w legal with respect to if there is a power A: e N such that w is a 
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subword of some realisation of Cm (a)- Furthermore, let ^ be the set of all legal 
words of length i with respect to i^,,,. We refer to the function C : N — > N, £ | i^^,, ^ | 
as the complexity function of C„,. It is known that the discrete hull of each member 
of Jfm has linear complexity, which implies that the topological entropy vanishes 
here. In the stochastic setting, the picture changes; see ID for background. 
Let m € N be arbitrary but fixed. The sets ;= 0, := {b}, := {a} and 

m in 

/=o;=o 

with denoting the Kronecker symbol, are called the generation sets of i^,,,. The 
product in ^ is meant to be the set-theoretic product with respect to concatenation 
of words. Moreover, we define := lim„^oo and refer to as the n-th generation 
set. The length of words in is given by the sequence €q = 0, = 1, ^2 = 1 and 
^„+i = ml„ +^„_i, for n>2. The set ^„ consists of all possible exact realisations 
Q~^{b). Since not all legal words result from an exact substitution, which can again 
be seen from the example bb, it is clear that |^^,| < C{£„) for n>2. 

In m, Godreche and Luck computed the topological entropy of under the 
implicit assumption that 

limllog(C(£„)) = limllog(|^^„|), 

which was recently proved by J. Nilsson IIT2I . This asymptotic identity is crucial 
because the exact computation of the complexity function of i^,,, is still an open 
problem. It is easy to compute | ^ | for £<m + 2 and it is known fTOl that 

I ^m^e I = t ( ■) - ^'«('« + 1 ) (3^ - 2m - 4) 

if m + 3 < £ < 2m + 2, while an extension to arbitrary word lengths seems difficult. 
In m, the entropy per letter for m = 1 is computed to be 

/ . li^ i^lM = E ^ - 0-444399 > 0, 
whereas a convenient representation for arbitrary m reads 

The result computed by Godreche and Luck for m = 1 can be recovered by the 
observation that (Aj — 1)/(1 — Af) = l/X^ in this case. Some numerical values are 
given in Table [1] It is not difficult to prove 1 10 1 that lim,„_!.oo /i,„ — 0, which can be 
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1 2 


3 


4 


5 


6 


7 




0.444399 0.408549 


0.371399 


0.338619 


0.310804 


0.287298 


0.267301 



Table 1 Numerical values of the topological entropy for RNMS with 1 < m < 7. 



verified by estimating the logarithm in (O via the square root and using the fact that 
A,,„/m tends to 1 as m 



4 Frequencies of subwords 



We adopt the method of computing the frequencies of subwords via induced sub- 
stitutions on words of length £ (with £ E N), which was introduced in ifTSl Section 
5.4.1], and modify it to fit the stochastic setting. To this end, we again assume that 
all probabilities in the definition of are strictly positive. 

If w = WflVVj • • • w^_i is a word of length £, we define w,- i to be the subword 



• • • Wj of w of length j ~ i+l. For £ > 2, we denote 
induced substitution defined by 



m,{i) ■ 



as the 



,(M) 



M+f-liyO<.<|C„,{H.('')|-l 



with probability Pj j , 



with probability p^ ^. , 



(3) 



where w*^'' € ( and v^'--'' is a realisation of w^'^ under with probability pj j. 
This way, we ensure that we are neither under- nor overcounting subwords of a 
given length. Similar to the case of i^,,,, the result is a random variable. 



The action of C 
and £^2: 



on words in f^^^ ^ is illustrated in the following table for m — 1 





w G 




y(''» P 


aa fl/jfl/j {ab){ba) p\ 


ab 


aba 


{ab){ba) p^ 


abba {ab){bb) p^Pi 




baa 


{ba){aa) Pq 


baab {ba){aa) p^Pi 


ba 


aab 


(aa) /?[ 


baba {ba){ab) p^ 




aba 


(ab) Pq 




bb 


aa 


[aa) 1 



Applying the lexicographic order to the words in Si^^^ f, leads to the corresponding 
substitution matrix M^^^ ^. For any fixed m € N and ^ = 2, we get 
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/{m-l)+pQp,„ (m-l)+Po 1-Pq 1\ 
^ ^-PoPm 1-/^0 Po 

"'.2 y-PoP,n 1 0- 

V PoPm 0/ 

This matrix has the spectrum (7(M,j^,) = {A„,, A/„, — /9q,/9qP„,}. Furthermore, it is 
interesting to observe that the spectrum of the matrices M,„ ( for £ > 3 is the same as 
that of M,jj 2^ except for the addition of zeros. Note that ^^^^ ^.j^ agrees with which 
implies that M,,, j = M. 

The substitution matrix M,,j ^ is primitive for all m and which allows an appli- 
cation of Perron-Frobenius theory; see (14] for general background. This implies 
that there is a strictly positive right eigenvector to the eigenvalue Note that 
A„j does not depend on any of the Pj at all, whereas this is not the case for 0^^'. 

We define a measure on the discrete hull X,„ as follows. For any word v S 
and ken, let Z^(v) := {w G X,„ | 

^[kk+e-i] ^ ^} '■^^ cylinder set of v that 
starts at position k. Then, the family {Z^,(v)}|(.gn generates the product topology and 
we define the measure fi : X„, R>q on the cylinder sets as /x(Z^,(v)) — ^'•'^'(v), 
where ^'^'(v) is the entry of 0^^' with respect to v. This is a proper (and consistent) 
definition of a measure on X,,,, which can also be found in ifTSl Section 5.4.2]. By 
construction, the measure is shift-invariant. 

The following theorem ITOl shows that, similar to the deterministic setting ifTSl . 
it is possible to interpret the entries of as the frequencies of legal subwords with 
respect to i^,„ as follows: 

Theorem 1 Let X,„ C be the hull of the random noble means substitution for 
m G N and jX the shift-invariant probability measure on X,„ as defined above. For 
any f & (X,„ , jl ) and for an arbitrary but fixed i e Z, the identity 

1 N+S-l n 

holds for jX-almost every x G X,„. 

The proof can be accomplished by inspecting the family of random variables 
R = {/(5'w)},gf^, where / is a patch recognition function that evaluates to 1 if 
{S'w)\^s,s+t~i] = ^ for an arbitrary but fixed word v G ^ and to otherwise. By 
observing that, given any ; e Z, the sets 

consist of pairwise independent words, we can split up the summation over R appro- 
priately and apply Etemadi's version of the strong law of large numbers Q Theorem 
1] to each sum over separately. This, in conjunction with an application of the 
Stone-Weierstrass theorem, implies the assertion. 
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5 Diffraction measure 

The symbolic situation is turned into a geometric one as follows. In view of the left 
PF eigenvector of M, a and b are turned into intervals of lengths A,„ and 1, respec- 
tively. The left end points are the coordinates we use. The corresponding continuous 
hull ¥„, is the orbit closure of all accumulation points of the geometric inflation rule 
under M. Let A C be a coordinatisation of an element of in M. Then, A is a 
discrete point set that fits into the same cut and project scheme as all elements of the 
family With respect to A, the smallest interval that covers A' = {x' | x g A} in 
internal space is given by [A/„ —1,1— A,',] . Then, A is relatively dense with covering 
radius A„,, and a subset of a model set, which implies that A is a Meyer set by ifTTl 
Theorem 9.1]. Let A^ = A n and consider the autocorrelation 

y := lim with y 8,. 

«^~vol(Bj « ^ 

The Umit almost surely exists due to the ergodicity of our system. By construction, 7 
is a positive definite measure which implies that its Fourier transform exists and is a 
positive measure. Regarding the Lebesgue decomposition 7= (7)pp + {j)ac + (7)sc' 
it is possible to compute the pure point part to be 

(7)pp= L \ria{-k') + U-k')\'8,, 

where L® = Z[A„,]/vn?+4 is the Fourier module. In the case of m— 1, the invari- 
ant measures 77^, 77^^ can be approximated via the recursion relation 

with « e N and ^ ;= Aj'. As ^" — > for n — > 00, an appropriate choice of the eigen- 
vector (na(0), 17^(0))^ for the equation 

V oyU(ov~ 'AW' 

which results from (01) for k ~ and n = 1, fixes the base of the recursion and 
provides the desired approximation. Since TJa(O) + J7z,(0) must be the density of 
A, which always is Aj/V5, one finds rjaiO) ^ l/VS and 77^,(0) = (A; - 1)/V5. 
Let /ij be the measure on induced by via suspension; see [S] Chapter 11] 
for general background. One consequence, due to a theorem of Strungaru IfTSl and 
an application of the methods of [SJ, is that our random dynamical system D := 
(Y„j,R,/i^) is ergodic, but not weakly mixing. In particular, it has strong long-range 
order. 
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Due to the stochastic setting with positive entropy, one expects a non-trivial abso- 
lutely continuous part. For m=l, the ergodicity of D almost surely yields a diffrac- 
tion measure which can be represented as 7 = (7)pp + (x{k) ■ X, where 

7(W) - liml.|E(g„(^))p. 

Here, E refers to averaging with respect to fi^ and 

aik) Jiml. (E(|g„(^)|V|lE(^„W)|') , 

with the random exponential sum g„{k) :~ LJiY e^^'^'*^^' and L,, := AjF,, + f'„_j, 
where is the n-th Fibonacci number. Now let 

A„{k) := E(^„(^)) and B„{k) E(|^„(^)|') - \E{g„{k)) \\ 
Godreche and Luck (S) derived a recursion relation for the sequence A„ (k), 

where Aq(A:) = ^-^'^^k Ay{k) — e^^Tni-^i Analogously, one derives a recursion 
relation for the sequence B^^ik), 

B„ik) = B„_iik)+B„_,ik)+2poPi •4„W, 

with 

A„{k) = (1 -cos(27r^L„_i)) ' |A„_2(^)|'+ (l -cos(2;r^L„_2)) • |A„_i 
-Re[(l -e^-^V.) . (1 _e-2-'^V2) .A„_,{k) 

and Bo(^) ^ ^i(^) ~ 0. In [81, almost surely by way of a misprint, the authors 
applied complex conjugation on A^^ ^{k) instead of A^^ which makes a huge 
difference, as the sequence B„(^) does not converge in that case. The recursion for 
B^{k) can be solved and the explicit representation reads 

u 

1=2 

A detailed discussion of the continuous part of 7 can be found in ifTOll . 

The illustration of an approximation of the diffraction measure 7 in case of m = 1 
and = /7[ = J, based on the sequences A„(A:) and B^{k), is shown in Figure [T] 
which agrees with the average over many realisations for the same length. 
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